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We investigate the finite temperature properties of the one-dimensional two-component Bose 
gas (2CBG) with repulsive contact interaction in a harmonic trap. Making use of a new lattice 
embedding for the 2CBG and the quantum transfer matrix we derive a system of two nonlinear 
integral equations characterizing the thermodynamics of the uniform system for all values of the 
relevant parameters: temperature, strength of the interaction, chemical potential and magnetic 
field. This system allows for an easy numerical implementation in stark contrast with the infinite 
number of equations obtained by employing the thermodynamic Bethe ansatz. We use this exact 
solution coupled with the local density approximation to compute the density profiles and local 
density correlation function of the inhomogeneous gas for a wide range of coupling strengths and 
temperatures. Our results show that the polarization in the center of the trap influences heavily the 
local correlator especially in the experimentally accessible Tonks-Girardeau regime. 
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I. INTRODUCTION 


As a result of the growing expertise in the creation and the manipulation of ultracold atomic gases we have witnessed 
in recent years the experimental realization of various strongly interacting many-body systems characterized by an 
extreme degree of purity and excellent control of relevant parameters including temperature, strength of the interaction 
and even dimensionality m- In particular, one-dimensional (ID) gases which can be created using optical lattices 
or atom chips are extremely interesting for several reasons. One reason is that these systems exhibit regimes such as 
the Luttinger liquid (LL) [4], spin-incoherent LL m\ and the ferromagnetic liquid [im, which are not present in 
2D and 3D. In addition, the high degree of control of the relevant parameters means that some of these systems can 
be well approximated by integrable systems providing a parameter-free comparison of theoretical predictions with 
measurements. The paradigmatic example in this case is the experimental realization [121117j of the Lieb-Liniger 
model [18] which has been investigated theoretically for more than fifty years. 

It is natural to expect that one-dimensional physical systems comprised of particles with internal degrees of freedom 
present a wider range of phenomena than their single component counterparts. Such spinor gases have been produced 
by trapping atoms in two or more internal states which can be referred as (pseudo)spin states. In the case 
of fermionic spinor gases the groundstate is antiferromagnetic [26j and at low-energy they are described by the LL 
theory. In addition the phenomenon of spin-charge separation is present which means that the Hilbert space separates 
in two independent sectors one for the collective spin excitations and one for the collective charge excitations. The 
situation involving bosonic spinor gases is more complex. While there are cases in which the physics is similar to 
their fermionic counterpart m, in the case of spin-independent interactions, which is integrable and will be the main 
subject of this paper, the groundstate is ferromagnetic [281129] . The spin excitations have a dispersion behaving like 
which makes the application of the LL theory impossible. For these systems the low-energy sector is described by 
a new universality class called the ferromagnetic liquid [^llj. 


In this article we are going to investigate the thermodynamics, density profiles and local density-density correlation 
function of the trapped ID two-component Bose gas interacting via a (^-function potential using the exact solution 
of the uniform system and the local density approximation (LDA). From the historical point of view, the first exact 
solution of an integrable multicomponent system was obtained by Yang m and Gaudin m for the two-component 
fermionic case using what we call nowadays the nested Bethe ansatz. The spectrum of the 2CBG was derived in [3^13^ 
and the low-lying excitations were investigated in |33| . The spin-wave excitations, spin dynamics, edge exponent in 
the dynamic spin structure factor, groundstate properties in a trap at T = 0 and impurity dynamics can be found in 
[33H111 • Even though the study of nonlocal correlation functions is still in its early stages [33 01] recent progress in 
the calculation of form-factors |43-|49| opens the way for the calculation of correlators along the lines of [13 El]. As 
a result of the integrability of the model one might believe that computing the thermodynamical properties should 
be an accessible task. Unfortunately, this assertion is not true. For example, the application of the thermodynamical 
Bethe ansatz (TBA) [52] [53] produces an infinite number of nonlinear integral equations (NLIEs) [54] which makes 
the extraction of physical information extremely difficult. Even though some important results can be derived on the 
basis of the TBA equations [5SHS7], it is highly desirable to obtain an alternative thermodynamic description allowing 
for an easy numerical implementation. Our result, Eqs. Q and Q, which was derived using the quantum transfer 
matrix (QTM) and a new lattice embedding of the 2CBG, provide such an alternative description. 

The basis of our method is the following observation. Gonsider an integrable continuum model and one of its lattice 
embeddings. By lattice embedding of a continuum model we understand an integrable model defined on a lattice 
whose Bethe equations and spectrum transform in a specific scaling limit in the Bethe equations and spectrum of the 
continuum model. Then, the thermodynamics of the continuum model can be derived from the thermodynamics of the 
lattice model if we perform the same scaling limit as in the case of the Bethe equations and spectrum. One example is 
the Yang-Yang thermodynamics [52j of the Lieb-Liniger model which was obtained from the thermodynamics of the 
critical XXZ spin chain [SH]- As we will show in Sect. IV in the case of the 2GBG the relevant lattice model is the 
critical q = 3 Perk-Schultz spin chain |59H64j . The thermodynamics of the spin chain is derived using the quantum 
transfer matrix method |65II69j . which can be defined only for lattice models but, has the fundamental advantage of 
providing a hnite number of integral equations. Performing the scaling limit in these equations we obtain Eqs. (§ 
and Q. This hnite system of equations coupled with the LDA allows for the calculation of the density prohles and 
the local density-density correlation function of the system subjected to a slowly varying harmonic potential. 

The plan of the paper is as follows. In Section |ll] we introduce the 2GBG and the thermodynamic description of 
the uniform system. The results for the density prohles and local correlator are presented in Section [ITT] Section [TV] 
contains the lattice embedding of the 2GBG and in Section |V] we remind some notions of the algebraic Bethe ansatz 
and introduce the quantum transfer matrix. The derivation of the free energy of the spin chain can be found in Sect. 
jVIj and the continuum limit is performed in Section |VII| The nested Bethe ansatz solution of the generalized g = 3 
Perk-Schultz model and the proof of several useful identities can be found in Appendices [^ and [^ Some of the results 
presented in this article were announced in the short paper [70]. 
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II. THE ONE-DIMENSIONAL TWO-COMPONENT BOSE GAS 


We consider a one-dimensional system of equal mass bosons with two internal degrees of freedom interacting via 
a (5-function potential and constrained on a ring of circumference Lb (periodic boundary conditions). The second 
quantized Hamiltonian in the presence of a harmonic trapping potential V{x) = miox^ is 


U 


Bose — 



2m 




I : : +(H(a;) - , 


( 1 ) 


with vl/ = ^ ^ 2 (^j ^ ^ 1 ^ 2 ( 2 ^) ) 5'2'z = ^ Q ^ ’ '"^here 'E'a(a;), a = {1,2} are ID quantum fields 

satisfying canonical commutation relations 'I>a(a:)'!'],(?/) — 'I'|,(?/)'Pa(a:) = SabS{x — y). In Q 5 > 0 is the coupling 
constant, /i is the chemical potential, H the external magnetic field (the Bohr magneton and the Lande factors are 
absorbed into H), uj is the trap oscillation frequency, m is the mass of the particles and : : denotes normal ordering. 
In experiments, the internal degrees of freedom are two distinguishable hyperfine states which can be thought of 
as a (pseudo)spin I. Compared with the scalar case (Lieb-Liniger model [TH]) the wavefunctions of the 2CBG are 
symmetric only under exchange of coordinates of particles with the same spin. One dimensional systems that are well 
approximated by the Hamiltonian 0 can be experimentally achieved in highly elongated cylindrical traps (w ^ a;j_) 
where uj± is the trap oscillation frequency in the transversal plane. Assuming that the 3D scattering length is 
much smaller than the transverse harmonic oscillator length l± = yjfi/muj± the coupling strength can be expressed 
as g ~ 2 / 1 ^ 031 )/(to/{^) = 2huj±a3D [7ll[72]. The realization of the ID regime requires that hu}± is larger than the 
thermal energy ksT and the chemical potential n [71 [71]. The experimental advances of the last decade paved the 
way for the experimental realization of such quasi-ID systems m characterized by coupling strengths which range 
from weak coupling (7 <C 1) to the strongly interacting Tonks-Girardeau regime (7 3 > 1). In order to lighten the 
notation in the following we are going to consider h = ks = 2m = 1 and introduce the effective interaction parameter 
c = gm/fi^. 


A. Bethe ansatz solution for the uniform system 


Even though the Hamiltonian Q is integrable only when the system is homogeneous iV{x) = 0), the trapped 
system can be efficiently investigated using the solution of the uniform system coupled with the local density approx¬ 
imation. Therefore, we will first review the exact solution of the homogeneous 2GBG and defer the treatment of the 


inhomogeneous system to Section HI Since we are interested in investigating the thermodynamic behavior we will 
consider only the case of repulsive interaction c > 0 (for c < 0 the system is thermodynamically unstable). For a 
system of M particles of which Mf are of type 1 and are of type 2 {M = -|- M^) the energy spectrum of 

the 2GBG obtained using the nested Bethe ansatz [351 [33] (see also [TSHTTj l is 


M 


Esose = ^eo{kf') - H{M^ - M 2 ), eo(fc) = - /i, 

i=i 


( 2 ) 


with satisfying the Bethe ansatz equations (BAEs) 
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(3b) 


The solution is characterized by two sets of rapidities, {kj^^},{k^^^} (a general characteristic of integrable two- 
component systems) with the second set of rapidities contributing to the energy 0 only via the BAEs. 

Even in the case of integrable models computing the thermodynamics is an extremely challenging task. One 
way to tackle this problem is the utilization of the thermodynamic Bethe ansatz |52[ 153] . In this framework, the 
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grandcanonical potential of the system is (/? = 1/T) |5^ 

f dk \n{l + T]i{k)), (4) 

J — OO 

with rji(k) satisfying the following infinite system of nonlinear integral equations 

ln77i(fc) = -/3(/c^ - fj.-H)+as ln(l + r]i{k)) + f * ln(l + ? 72 (fc)), (5a) 

ln77„(fc) = / * (ln(l + ry„_i(/c)) + ln(l + ry„+i(fc))) , n = 2, • • • , oo , (5b) 

together with the asymptotic condition lim„_).oo ln? 7 „(/c)/n = 2(311. In Eqs. (j^ g * h(k) = g(k — k')h(k') dk', 
f{k) = l/[2ccosh(7rfc/c)] and a„(/c) = nc/[27r((nc/2)^ + A:^)]. While important, it is obvious that extracting physically 
relevant information from this system of equations is very hard even from the numerical point of view (see [551457) 1. 
highlighting the need for a more manageable thermodynamic description of the 2CBG which will be presented in the 
next section. 


B. Alternative thermodynamic description of the uniform 2CBG 


A more efhcient thermodynamic description of the 2CBG was proposed in IZO]. The full derivation of this result 

which is based on the connection of our model with the (-) Perk-Schultz spin chain and the quantum transfer 

matrix method represents one of the main results of this paper and will be presented in Sections VI and VII In this 
description the grandcanonical potential per unit length of the 2GBG is 


^{d-^H,j 3 ) = [ [ln(l + ai(A:))+ln(l + a2(A:))] dk, 

Jr 

with ai, 2 (^) auxiliary functions satisfying the following system of integral equations 

\naiik) = -/ 3 {k^-g-H)+ [ (k - k')\n{l + ai{k')) dk'+ f K^ {k - k'- ie) \n{l + a2ik')) dk', 

Jr Jr 

\na2ik) = -l 3 {k^ - H + H)+ [ Kf{k-k' + ie)\n{l + aiik'))dk'+ f {k - k')\n{l + a2{k')) dk', 

Jr Jr 


where £ —)• 0 and 


Kl^ik) = 


^ K^(h') = ^ ^ K^(k) = ^ 

27 t k"^ + ^ 27r k{k + ic) ' ^ 27r k{k — ic) 


( 6 ) 


(7a) 

(7b) 


( 8 ) 


We can check the validity of our result in three well known limits. First, we will address the noninteracting limit, 
c —>■ 0. Using limc^o K^ik — ie) = limc_>.o Ki{k + ie) = 0 and limc_j.o Koi^ ~ ~ NLIEs Q decouple 

with the result 

Inai(fc) = —I3{k‘^ — d- — H) + lii(l + ai(fc)), 
lna 2 (fc) = —P{k'^ — g, + H) + ln(l + a 2 {k)). 

These equations are easily solved obtaining 


<^(m, H,(3) = J- I ln(l - + in(I - 

Jr 1 


dk , 


which is exactly the grandcanonical potential of two noninteracting Bose gases at different chemical potentials. 

In the strong magnetic field limit, H oo, due to the fact that the inhomogeneity —/3(A:^ ~ P + H) is large and 
negative 02 (fc) ~ 0 , which means that the thermodynamics of the system is given by 

(l){g,H,(3) = -J— [ \n{l + ai{k)) dk , \ogai{k) =-j3{k^ - g - H) + [ {k - k')\n{l + ai{k')) dk'. 

2'7rp Jr Jr 


^ This form of the TBA equations can be easily derived from the one used in [SIET] by the use of simple transformations and identities 
as in Chap. 12 of [53]). 
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FIG. 1. (Color online) Plot of the grandcanonical potential per length of the 2CBG and impenetrable 2CFG (dot-dashed line) 
as a function of temperature for /j, = 1 ,H — 0.5 (left panel) and fj, = 100 , H = 80 (right panel). For large values of the coupling 
constant the 2CBG result approaches asymptotically Eg. ([^ . Grandcanonical potential per length, temperature, chemical 
potential and magnetic field in units of t/io, To, Mo and ho I78| . 


This result, which is the Yang-Yang thermodynamics [55] of the Lieb-Liniger gas, confirms the natural expectation 
that in the strong magnetic limit the system will become fully polarized, and, therefore, the thermodynamics will be 
the same as in the single component case. In a similar fashion, we can show that if the magnetic field is fixed and 
finite, then, in the low temperature limit (T ^ /x, c), the same formulas are obtained proving that the ground state 

is ferromagnetic (fully polarized). 

In the limit of impenetrable particles, c —>■ oo, the system is effectively “fermionized” and we should reproduce the 
result obtained by Takahashi (Chap. 12 of [53]) for impenetrable repulsive two-component fermions (2CFG) 

= [ *ln(l + 2cosh(/3iT)e-'5('='-^)) . (9) 

Even though we did not succeed in providing an analytic proof of this equivalence, we have checked numerically and 
found perfect agreement as can be seen in Fig.[^ 


C. The “heuristic” derivation 

The rigorous and self-contained derivation of the thermodynamic description presented in the previous section will 
be given in Sections | VI| and [Vll| Here we provide a heuristic derivation based on the fact that the system of equations 
^ share the same structure with the NLIEs for the spin ^ XXX spin chain. Let us be more precise. We consider a 
very general XXX spin chain case characterized by the Hamiltonian Hxxx = TTo + PnJn with nearest-neighbor 

interaction Hq and higher conserved currents l7n- Explicitly we have Ho := »+i + h)-^T,i S! , with J 

and H the interaction strength and magnetic held. The TEA result for the free energy per unit length is E(iL,/3) = 
eo - dkf(k) ln(l + fji(k)) , (cq is the zero point energy and f(k) = l/[2 cosh(7rfc)], / = / + ^i^h 

7)1 (fc) satisfying 

Inrji(k) = -2xj3jj{k) + f * ln(l -b mik)) , (10a) 

ln7)„(A:) = / * [ln(l-b ? 7 n_i(/c)) -b ln(l-b ? 7 „+i(fc))] , n = 2 ,---,oo, (10b) 

together with the asymptotic condition lim„_>.oo ^'Rfin{k)/n = j3H. A more compact result can be obtained using the 
quantum transfer matrix [55] [Ml ES] yielding F{H,I3) = Cq — dk f{k)[\n{l + di(fc)) -b ln(l -b a 2 (^))], with 
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only two auxiliary functions satisfying 

Inai(fc) = —/3[27rJ/(fc) — H/2] -\- Kq * ln(l + Oi(fc)) — K 2 * ln(l + a 2 {k )), (Ha) 

lna 2 (fc) = —/3[27rJ/(fc) + H/2] — Ki * ln(l + oi(fc)) + Kq * ln(l + a 2 {k)). (Hb) 


where 


Km = ^ 


r+00 ^ikx^-\x\l 1 


„ , ^ dx, Ki(k) = — 

2 cosh § ’ ^ ^ 27r 


r*+oo ^ikx^x—\x\/2 


1 


dx, K 2 {k) = — , 
2 cosh I 27r J 


+00 ^ikx^—x—\x\/2 


2 cosh ^ 


dx. 


Comparing the two systems of equations (10) and (11) and their corresponding expressions for the free energy we 
notice that a) ln(l + r 7 i(fc)) = ln(l + ai(fc)) + ln(l + a 2 (^) and b) the driving terms (modulo the magnetic field) in the 
QTM system are the same as the driving term of the integral equation for \nfji{k). Due to the similar form of the TEA 
equations for the XXX spin chain (10) - with practically general function / - and the 2CBG ([^ it is natural that a 
similar system like 0 can be derived via: a) fc —>• fc/c b) ln(l+77i(fc)) = ln(l+ai(/c))+ln(l+a 2 (^) and c) replacement 
of the driving terms with — /3{k‘^ — H) + as* f *\n{l + rii{k)) = ——/r=Fi7) + a 3 */*(ln(l + ai(fc)) + ln(l + a 2 (A:))). 
Performing these transformations we obtain (|^ for the grandcanonical potential and the following system of integral 
equations 


\nai{k) =-I3{k'^ - fi +H) + {as* f + Kq^ , 

lna 2 (/c) = —/?(fc^ — fi — H) + {as * f — Ki) * ln(l + 


m, 


,) * ln(l + di{k)) + {as* f - K 2 ) * ln(l + 02 

oi(fc)) + (as * / + Xo) * ln(l + a2(fc)) 


( 12 a) 

( 12 b) 


Finally, it is easy to show that {as*f + Ko){k) = KqW ) (a 3 */~^ 2 )(^) = Ki{k + ie), {as* f — Ki){k) = Ar^(A: —je) 
by taking the Fourier transform and applying the convolution theorem, completing the derivation. 


III. DENSITY PROFILES AND LOCAL CORRELATION FUNCTIONS 

The experimentally relevant case of a trapped system {V{x) ^ 0 in 0 ) can be investigated using the exact solution 
of the uniform system, Eqs. ([^ and ([^, and the local density approximation. Under this approximation, which is 
valid for a slowly varying potential, the system in a trap can be described locally m as a uniform gas with chemical 
potential and magnetic field defined by 


^i{x) = fi{0) - V{x), H{x) = H{0), 


(13) 


with ^(0) and H{0) the chemical potential and magnetic field in the center of the trap. Therefore, for given values of 
temperature T, coupling strength c, chemical potential and magnetic field in the center of the trap /i(0) and H{0), 
relevant thermodynamic quantities at a distance x from the center of the trap can be computed using Eqs. ([^ and 
0 with ^ and H replaced by ( [T^ . 

In this section we will be mainly concerned with the calculation of the density profiles in the trap and the normalized 

('~p\ 

finite temperature local density-density correlator ' [731 EDI [S]. For a uniform system in thermal equilibrium the 
linear densities of the two types of bosons can be obtained from the derivatives of the grandcanonical potential per 
unit length 0 : 


ni = 


1 ( d(j) d(j) 


2 \dg dH 


n2 = 


d(f) 

dg 


d^\ 

dH) ■ 


(14) 


In the case of the local correlator, as it was shown in a simple application of the Hellmann-Feynman theorem 
allows to express this important observable as the derivative of the grandcanonical potential with respect to the 
coupling strength 


(T) 

92 


{x) = 


J2a,bi'^iix)'^lix)'^aix)'i>b{x))T 

{Ea{'^U^)'^a{x))T 


d(j) 


(Ea’^o)^ dc ■ 


(15) 


The local density correlation function of the uniform system which presents a nonmonotonic behaviour as a function 
of temperature was investigated by Caux, Klauser and van den Brink in |561157j . 

A uniform 2CBG can be described with the help of three dimensionless parameters: interaction parameter 7 = 
c/{ni -I- 712 ) = c/n, reduced temperature r = T/T^ = T/n? and polarization P = {ni — n 2 )){ni + 77 , 2 ). While this 
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FIG. 2. (Color online) Diagram of the different regimes of the fully polarized uniform 2CBG (Lieb-Liniger model) in the 
7 — t plane (boundaries depicted as dashed lines) [73]. The six regimes depicted are defined as follows: Tonks-Girardeau 
(TG) 7 S> 1 ,t 7 “^; High-Temperature Tonks-Girardeau (TG-high T) 7 S> 1,7~^ <C t <C 1; Gross-Pitaevskii a (GPa) 
7 ^ l,t <C 7 “^; Gross-Pitaevskii b (GPb) 7 “^ <C t <C Decoherent Quantum (DQ) 7 “®/^ <C t ^ 7 “^; Decoherent 

Classical t ^ max{'y~‘^, 1). The black dots are the values of 7 ( 0 ) and t in the center of the trap of the representative systems 
considered in the next section with the continuous black lines at t = {5 x 10“^, 5 x 10“^, 1, 5 x 10^, 5 x 10^} representing the 
values of 7 ( 0 :) in the trap for these systems. 


description can be naturally extended to the inhomogeneous case by replacing the density n with the local value in 
the trap n{x) it is preferable to work (see m) with 


7(x) 


c 

n{x) ’ 



P{x) 


ni{x) — n 2 {x) 
ni{x) + n2{x) 


(16) 


where we have introduced a new temperature parameter t = t{x)/^‘^{x) which has the advantage of not depending 
on the local density n{z) and characterizes the entire system in thermal equilibrium. 


A. Regimes in a trapped and uniform 2CBG 

Analytical results on the local density correlation function of the 2CBG at arbitrary polarization (any value of H) 
are very scarce in the literature (a notable exception is EZ])- However in the case when the system is fully polarized 
(P = 1) the 2CBG is equivalent to the Lieb-Liniger model for which a classification of different regimes exist. The 
fully polarized case represents a good starting point for the investigation of the 7 — t — P parameter space if we 
take into account that ID bosons favour ferromagnetic behavior for spin independent interactions. In the case of 
the Lieb-Liniger gas, Kheruntsyan, Gangardt, Drummond and Shlyapnikov IZIISII identified three regimes (each 
containing two sub-regimes) using the properties of the local density correlator as a function of coupling strength and 
temperature. Using 7 and t as parameters this regimes can be characterized as follows: 

Strong coupling regime. In the limit of strong interaction (or low-densities), 7 ^ 1 , the system becomes equivalent 
with a system of free fermions and therefore the local density correlation function is suppressed, g 2 <SC 1 , as a 
consequence of the fermionic nature of the wavefunctions. Two sub-regimes are distinguished: a) Tonks-Girardeau 
characterized by 7 ^ 1 and temperature smaller than the degeneracy temperature t ^ 7 “^ and b) High-Temperature 
Tonks-Girardeau for which 7 “^ <C t <C 1. 

Gross-Pitaevskii regime. In the limit of vanishing interactions (or large densities), 7 <C 1, the local density corre- 
lation function approaches g^ — \ which is the value for free bosons. The two subregimes are: Gross-Pitaevskii (a) 
(7 <C 1 7 “^) and Gross-Pitaevskii (b) ( 7 “^ <^t 7 “^/^). 
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7(0)=1.7x 10“^ t=500 


7(0)=9.7x 10“^ t=500 


7(0)=2.2x 10“^ t=500 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 0 

1.0 




X/^GP 


X//?GP 


FIG. 3. (Color online) Density profiles (upper panels) and local density correlation function g^\x) (lower panels) of a trapped 
2CBG at temperature t = 5 x 10^ for three values of 7 ( 0 ) = {1.7 x 10“®,9.7 x 10“^, 2.2 x 10“^} and various polarizations. In 
the upper panels the thick continuous, thick dashed, and thick dot-dashed (thin continuous, thin dashed and thin dot-dashed) 
lines represent the majority (minority) components and n(0) is the total density in the center of the trap. The insets depict 
the variation of the polarization in the trap. 


('T\ 

Decoherent regime. At high temperatures the local pair correlation is close to 52 ' — 2- The two subregimes are: 
Decoherent Quantum ( 7 “^/^ ^ t <C 7 “^) and Decoherent Classical {t :$> max{^~'^, 1)). 

For a trapped system this classification still remains valid with the possibility of different regimes coexisting within 
the trap. Due to the fact that for large distances from the center of the trap the density vanishes n{x) —>■ 0, 'y{x) 00 

this means that the tails of the system are always in the TG regime (for t < 1) or the DC regime (for t > 1). For given 
values of the interaction and temperature parameters 7 ( 0 ) and t the density profile follows a straight line at t parallel 
with the 7 axis as it is shown in Fig. [^for t = {5x 10“^, 5 x 10“^, 1, 5 x 10^, 5 x 10^}. All these considerations are valid 
for the fully polarized 2CBG. In the following sections we are going to investigate the influence of the polarization on 
the density profiles and local density correlation function. 


B. Density profiles and local correlation functions at high-temperatures 


For a given value of the coupling strength c the region of high temperatures is defined by f 3> 1. We have calculated 
the density profiles and the local correlation function on a fine grid of the trap using Eqs. ([^, Q, ( [l4| ), (151 and 
the LDA chemical potential and magnetic field (13). The results for three values of 7 ( 0 ) at f = 5 x lO^and various 
polarizations are shown in Fig. In Fig. the distance x is plotted in units of the Thomas-Fermi radius in the Gross- 
Pitaevskii regime ca, Rgp = ( 6 iVc/a;^)P^ with N the number of particles in the trap N = f(ni(x) + n2(x))dx, 
and uj is the trap oscillation frequency. The left panels of Fig. depicts the situation in which the gas in the center 
of the trap is in the GPa regime with three polarizations very close to 1 (see Fig. where 7 ( 0 ) = 1.7 x 10“^ is the 
leftmost point on the t = 5 x 10^ line) for which ~ 1. While in the center of the trap the local correlator 

is almost constant showing no dependence on P(0) the situation changes dramatically in the tails where the gas is 
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7(0)=7x10'^ 1=5x10“ 



0 2 4 6 8 10 12 

X//?GP 


FIG. 4. (Color online) Density profiles (npper panel) and local density correlation function g^\x) (lower panel) of a trapped 
2CBG at temperature t = 5 x 10“ for 7 ( 0 ) = 7 x 10“® and polarizations P(0) = {0.0003,0.6918,0.9999}. In the upper panels 
the thick continuous, thick dashed, and thick dot-dashed (thin continuous, thin dashed and thin dot-dashed) lines represent 
the majority (minority) components and n(0) is the total density in the center of the trap. The inset depicts the variation of 
the polarization in the trap. In the lower panel the thin lines denote the asymptotic values given by Eq. (|17[) 


in the DC regime. Here we see that even small variations of the center polarization have a significant effect in the 
tail behavior of the local correlator and polarization. For larger values of 7 ( 0 ) (center and right panels of Fig. the 
effect of the polarization in the center of the trap on the local correlation function becomes more pronounced with 
g^\x) increasing both as a function of 7 and P{0). In addition we see that the density profiles present larger tails 
as we approach the DC regime. For large values of 7 and t, when the 2CBG is in the DC regime, it can be shown 
using Wick’s theorem and Boltzmann distribution [sni ET] that for a uniform system the asymptotic values of the 
local correlator and polarization are 

= _Z_r_ _ 7 _ (■17') 

^2 ^ g/3(/x+i/) g/3(A*-ff) • GC 

This means that in the case of a trapped system at very large temperatures the local correlator at the edge of the 
sample should tend asymptotically to 2 for P(0) = 1 (very large magnetic field) and to 1.5 for P(0) = 0 (vanishing 
magnetic field). This can be seen clearly in Fig. where we present the density profiles and local correlator for 
t = 5 X 10"^ and 7 ( 0 ) = 7 x 10“^. Here the entire sample is in the DC regime and the effect of P(0) on is very 
large. 


C. Density profiles and local correlation functions at intermediate temperatures 

For temperatures close to 1 we can encounter the situation in which for small variations of t the profile of the local 
correlation function changes drastically as it can be seen in the left panels of Fig. for the fully polarized gas with 
7 ( 0 ) = 0.11 and t = {1, 2, 3}. While the density profiles are almost identical, g^ {x) presents a highly nonmonotonic 
behavior with asymptotic values at the edge of the samples varying from 0.65 for < = 1 to 1.05 for f = 3. A similar 
phenomenon can be seen for fixed temperature and various polarizations in the center and right panels of Fig. For 
t = 1 and 7 ( 0 ) = 0.5 the local correlation functions are almost equal in the center of the system but for large distances 
g 2 \x) is suppressed for small polarizations and enhanced for large polarizations. The role of P(0) increases at strong 
coupling (large 7 ( 0 )) as shown in the right panel where for 7 ( 0 ) = 1.94 and t = 1 the local correlators present a strong 
dependence on the polarization even in the center of the trap. 
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7(0)=0.11 t=1,2,3 


7(0)=0.50 t=1 


7(0=1.94 t=1 


0.0 0.2 0.4 0.6 0.8 1.0 1.2 0.0 0.5 1.0 1.5 0.0 0.5 1.0 1.5 2.0 



X//?GP 


x/Rgp 


X//?GP 


FIG. 5. (Color online) Density profiles (upper panels) and local density correlation function g^\x) (lower panels) of a trapped 
2CBG at different temperatures and polarizations for three values of 7 ( 0 ) = {0.11,0.5,1.94} and various polarizations. In the 
upper panels the thick continuous, thick dashed, and thick dot-dashed (thin continuous, thin dashed and thin dot-dashed) lines 
represent the majority (minority) components and n(0) is the total density in the center of the trap. The insets depict the 
variation of the polarization in the trap. 


D. Density profiles and local correlation functions at low-temperatures 

At low-temperatures (t <C 1) the edges of the system will always be in the TG regime and, therefore, the local 
correlator will be strongly suppressed. In addition, we expect the density profiles to be sharper as it can be seen 
in Fig. for t = 5 x 10“^. As expected g^\x) decreases in the tails of the distribution until it reaches a limiting 
value which depends on 7 ( 0 ). The influence of the center polarization becomes more pronounced at stronger coupling 
with AO) monotonically increasing as a function of P(0). It should be noted that numerical investigations at 
low-temperatures and low-polarizations {H —)■ 0) are very difficult due to the fact that in this regime the NLIEs 0 
become numerically unstable (the same phenomenon happens in the case of the TEA equations®) as a consequence 
of the first-order phase transition at T = 0,= 0. For t = 5 x 10“'^ results are shown in Fig. [7® 

The numerical data presented above show that the polarization at the center of the trap has a profound influence 
on the density profiles and the local density correlation function of the inhomogeneous system. g^\o) is a monotonic 
increasing function of P(0) with the influence of the center polarization being more pronounced at strong coupling. 
Also, the local correlation function which takes a wide range of values between 0 and 2 serves as a better discriminant 
than the density profiles for the different regimes of the 2CBG. 


IV. THE 2CBG AS THE CONTINUUM LIMIT OF THE PERK-SCHULTZ SPIN CHAIN WITH (-) 

CRADINC 

Our method of deriving the thermodynamic description of the 2GBG, Eqs. § and 0 , is based on the fact that 
the continuum model can be obtained by performing a specific limit of an appropriate lattice model for which the 
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7(0)=0.61 t=5x10 


,-2 


r(0)=3.68 t=5x10 


-2 


y(0)=8.17 t=5x10 
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FIG. 6 . (Color online) Density profiles (upper panels) and local density correlation function g^\x) (lower panels) of a trapped 
2CBG at temperature t = 5 x 10“^ for three values of 7 ( 0 ) = {0.61, 3.68, 8.17} and various polarizations. In the upper panels 
the thick continuous, thick dashed, and thick dot-dashed (thin continuous, thin dashed and thin dot-dashed) lines represent 
the majority (minority) components and n(0) is the total density in the center of the trap. The insets depict the variation of 
the polarization in the trap. 


quantum transfer matrix technique can be employed. For the two-component Bose gas the relevant lattice model 
is the q — 3 Perk-Schultz spin chain. In this section we are going to present this scaling limit and show how the 
grandcanonical partition function of the 2CBG can be derived from the canonical partition function of the spin chain. 
The Hamiltonian of the q = 3 Perk-Schultz spin chain with arbitrary grading is |61H64) 


'H'PS = J£ 


L 

E 


3 

a—1 


E< 

a,b—l 

a^b 


-'ab ^ba 


I Sin 7 


3 

E 

i,b—l 

a^b 


sign(a - 


\ 


/ 


-EE^° 

j=la=l 


oU) 


(18) 


where L is the number of lattice sites, 7 G ( 0 , 7 r) determines the anisotropy, {q = e®'*'), J > 0 determines the strength 
of the interaction and ha are chemical potentials. The parameters { 61 , 82 , 63 ) can take the values ±1 and we will 


call them the grading of the system. In (18), e 

{6ab)ij — hf 


U) _ 

ab 


7-1 


^ab ' 




with Cab the 3-by-3 matrix with elements 


Sbj and I 3 the 3-by-3 unit matrix. The Perk-Schultz Hamiltonian is the sum of Hb (in the brackets) which 
is the fundamental spin-model associated with the trigonometrical Perk-Schultz R-matrix |62j and hLc the chemical 
potential part which does not break the integrability. In the next section and Appendix we will show that for 
(ei,e 2 ,e 3 ) = (-) the energy spectrum is 


M 


Eps = ^eo(u7) +Mi(/i 2 - hs) + Eq , eo{v) = J 
1=1 


. 2 

Sin 7 


sin(z; — 7 ) sin 7 ; 


+ — /i 2 , Eq = JL cos 7 — hiL . 


(19) 
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y(0)=7.37 t=5x10”'' 


y(0)=19.29 1=5x10““* 


y(0)=88.41 1=5x10““* 


0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.0 0.1 0.2 0.3 0.4 0.5 



FIG. 7. (Color online) Density profiles (upper panels) and local density correlation function g^\x) (lower panels) of a trapped 
2CBG at temperature t = 5 x 10“"^ for three values of 7(0) = {7.37,19.29, 88.41} and various polarizations. In the upper panels 
the thick continuous, thick dashed (thin continuous, thin dashed) lines represent the majority (minority) components and n(0) 
is the total density in the center of the trap. The insets depict the variation of the polarization in the trap. 


with satisfying the BAEs 


sin(7 — 

(1) 

sint!{ 


sin('(;i^^ - vf'‘ - 7 ) ^ ^ sin(?;i^^ - v^p'’ + 7 ) _ ^ 

11 . / (1) (T)—; 7 111“'* , s-i," 


f=i sin('(;}^'' - + 7 ) ' ' sin(r;i^* - 4 ^*) 


,M, 


“ sin(vp* - - 7 ) sin(t-p* - 4 "^ - 7 ) 

111 ) , , ( 2 ) ( 1 ), 11 __( 2 ) , 


1=1 


1 = 1 . 


. / {Z\ (i)\ AJ. . / IZ) 

sva[v\'-v)’) p=iSin(i;}'- 


,Mi. 


( 20 a) 


( 20 b) 


Up 


p^i 


First, we will show how the BAEs @ of the 2CBG can be obtained from the BAEs ( [20| ) of the Perk-Schultz spin 
chain. The spin chain is characterized by the following set of parameters: number of lattice sites L, anisotropy 7 , 
strength of the interaction J > 0, lattice constant <5 and chemical potentials hi,h 2 ,h^. Consider —>• / 1 ^/2 


and 


,( 2 ) 


iSkf'^/e + 7 + 7 r /2 with e —>■ 0. The BAEs (20l become 


/ sinh((5fci^*/e - 17 / 2 ) \ si-ah{5k^p/e - 5k^p/e - i^) cosh(i5fci^Ve - Sk^^^/e + i^/2) 

\ sinh((5fci^*/e + 17 / 2 ) / smh{Ski^'^/e — Sk^^^/e + i^) cosh{Ski^^/e — Sk^,^^/e — i^/2) 

1 /s 

cosh((5fcp*/e — + 17 / 2 ) ^ sinh(( 5 fcp'*/e — + * 7 ) 

cosh((5fcp*/e — — 17 / 2 ) sinh{6kl‘^^/e — Sk^^/e — i^) 

p^i 
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TABLE I. Parameters for the (-) Perk-Schultz spin chain and the spinor Bose gas. 


(-) Perk-Schultz spin chain 

Spinor Bose gas 

lattice constant S — >■ 0{e^) 
number of lattice sites L — > Cl(l/e^) 
interaction strength J > 0 
chemical potential hi — >■ 0 {e^) 
chemical potentials /12, ha — >■ 0 {t'^) 
inverse temperature /3 
anisotropy 'y — tv — t 

particle mass m = 1/2 
physical length Lb = LS 
repulsion strength c = /S 

chemical potential fi = -+ 4 ^ + ^(ha -I- ha) 

magnetic field H — (ha — h 2 )/{ 2 S^) 
inverse temperature /3 = 


In the second step we perform 7 = tt — e with the result 


/ cosh(<5A:i^Ve + fe/ 2 )\ sinh(Ski^^/e — Skj^^/e + ie) sinh(5fcs^^/e — — fe/ 2 ) 

\cosh((5A:i^V^ ~ *^/ 2 ) / j=i sinh(Ski^^/e — Sk^^/e — ie) sinh((5fcs^V£ ~ Sk^'^/e + ie/2) 

i/s 

sinh(i5A:| /e — Skj /e — *e/2) .pj- sinh(i5A:pV^ ~ Sk^^'^/e — it) 

j=i sinh(i5A:p^/e — 6kj^^/e + it/2) p^i sinh(i5A:p^/e — 5k^'’/t + it) 

p^i 


( 21 a) 


(21b) 


Finally, taking the limit L —i- 00 like 0(l/e^), 5 —>■ 0 like C>(e^), such that Lb = L5 and c = t^/5, equating Mi = Mi 
and using 


cosh{S ki^^/t + it/2) l + iki^^5/2 
cosh{5ki^'^/t —it/2) 1 — iki^^5/2 


we find that (21 1 transform into the BAEs ^ of the Bose gas. Performing the same transformations in the expression 
for the energy we obtain 


M 

Eps -Eo = Y^ V - - J^V^] + (hi - h2)M + {h2 - h3)Mi + 0(e®), 


i=i 


( 22 ) 


which shows that in order to obtain the energy spectrum of the Bose gas ^ we need to scale J —>■ 00 like C>(l/e^) 
and hi —>■ 00 like 0(l/e^) such that J6^ = 1 and Jt^ + Jt^/4: — hi + h 2 is finite. Therefore, in the thermodynamic 
limit performing this scaling we have 

e'5^'>Z(hi,h2,h3,/3)^Z(/i,iJ,/3), (23) 


where /? = l/T, Z(hi, h 2 ,h^, (3) is the canonical partition function of the spin chain and H, /?) is the grandcanon- 
ical partition function of the 2CBG. 

This scaling limit is independent of the temperature of the systems. However, taking into account that we are 
interested in the thermodynamic behavior it is more useful to introduce a scaling limit which involves also the 
temperature. This is justified because from the the rmo dynamics point of view the energy spectrum enters via the 


expression e . Multiplying with /3 the identity (22) we find that (3[Eps — Eq) = [3EBose, with [3 the inverse 


temperature in the continuum system, if /3 = /3/(5^, J = 1 hi —^ 0 like 0{t^), such that — hi/S"^ is finite and 


^ 2,^3 —>• 0 like 0(t^). In this case Eq. (23) becomes 


e^^°Z{hi,h2,h-i,l3) Z{^i,H,j3). 


(24) 


In Table □we present in a compact form this continuum limit (the spectral parameter in the lattice model scales like 
V —> iSk/t) which will be used in the next sections to derive the thermodynamic behavior of the 2CBG from similar 
result for the (-) Perk-Schultz spin chain. 

















V. ALGEBRAIC BETHE ANSATZ AND THE QUANTUM TRANSFER MATRIX OF THE (-) 

PERK-SCHULTZ SPIN CHAIN 
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In the next section we are going to derive a finite set of NLIEs characterizing the thermodynamic of the (-) 

Perk-Schultz spin chain using the associated quantum transfer matrix. Therefore, it will be useful to recall some basic 
notions of algebraic Bethe ansatz |82) . 


A. ABA for the Hamiltonian 


The Hamiltonian (18) represents the fundamental spin model (Chap. VI of 

T)H] defined by 


]) associated with the trigonometric 


q = 3 Perk-Schultz R-matrix 


R(w,ri;) = ^ R““(u, w)eaa d Caa + ^ w) Boo <8 e&b -I- ^ i w) Cab eba , 


(25) 


a,b—l 

a^b 


a,b—l 

a^b 


withE] 


sin [7 + eaiv - w)] 


R-(„, , Kliv, W) = ^ ( 26 ) 


smy 


a^b sm 7 


a^b 


and (eab)ij = ^ia^jb the canonical basis in the space of 3-by-3 matrices. This 9-by-9 R-matrix has the property that 
R(0, 0) = P, with = (^Qjbjdaabi the permutation matrix, and satisfies the Yang-Baxter equation 


a'^b',c' — l a',b',c' — l 


(27) 


Let us sketch how we can obtain the Hamiltonian (18 1 (more precisely, the component within the brackets) in the 
framework of ABA. First, we need to introduce the 


--operators defined by 
L,(u,0) = ^ Rn)(v0)ei°)e« ^ , L,(u,0) G End ((C3)®(^+i)) 


(28) 


a,b,ai ,61 = 1 


where is the canonical basis of operators acting on (C^)®('^+^\ i.e., = Cab 0 and = I 3 0 0 

Cob 0lf for j = 1, • • • , L. These operators act on the tensorial product of C^, which is called the auxiliary space, 
and the Hilbert space of the spin chain S) = (C^)®'^. In this auxiliary space, using the definition of the R-matrix (251, 
the L-operators can be represented as 


Lj(u,0) = 


7+(u,0)e^^J 
7 +(u, 0 )e 


'U) , Ai) 

^22 ^ ^33 


7-(v0)e^^ 


a2(^^,0)e^'’2^ -b/3(u,0) 


Sj) I So) 
^11 ^ ^33 


(i) 

13 


(i) 

23 


7_(u,0)e^^i^ 
0) 


\ 


where we have introduced 


7 +(w, 0 )e; 

sin(u — w) 


7_(u, 0)6^2 
a3(i',0)e^'^3^ -b/3(u,0) 


SA , AA 


(29) 


= sin [7 + g.(u-u;)] ^ ^ snq^^ ^ ^ ^±.,„(a-b)(.-.) ^ ( 30 ) 


Sin 7 


81117 


In (29), now represents the canonical basis of operators acting on and, therefore, the elements of the 


matrix represent operators acting on the Hilbert space of the spin chain. The monodromy matrix is the ordered 


^ In |621163| the authors considered the more general case = Gab 7 ^^ with GabG^^ = 1 (no summation). The situation 

considered in our paper corresponds to Gab = 1 




















product of L-operators 
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Tlliv) Ti2{v) Ti 3 {v) 

T(u) = LL{u,0)Li_i(u,0) • • • Li(u,0), ^(u) = | T2i(u) T22{v) T23(w) 

T3i{v) T32{v) 733 (w) 

and provides a representation of the Yang-Baxter algebra 

R(u, iu)[T(u) 0 T(u>)] = [T(u) 0 T(r(;)]R(u, w), 


( 31 ) 


(32) 


where Rj,^ j^(u,w) = (PR)j(^^(w,re) = (u, w). In Eq. (32| the tensor product should be understood as the tensor 


product of 3-by-3 matrices with operator valued entries as presented in the r.h.s of (31). Finally, the transfer matrix 
is defined as the trace of the monodromy matrix in the auxiliary space 


t(u) = troT(z;) = rii(u) + T22(u) + T^^iv). 


(33) 


Following [83] we can show that 


L 3 


T-lps = Jeisinjt ^(O)t'(O) — 


,U) 


(34) 


j=l a=l 


Obtaining the eigenvalues of the transfer matrix requires the existence of a pseudovacuum on which the monodromy 
matrix acts triangularly. Using (29), it is easy to see that 



(35) 


L times 


satisfies the requirements of a pseudovacuum, and 


T(u)|0) = 


/ rn(u)|0) Ti2(u)|0) Ti3(u)|0) 

T2liv)\n) T22iv)\n) T 23 iv)\n) I = 

\T3,iv)\n) T32(u)|0) T33iv)p 


( (ai(u,0))^|O) ri2(u)|0) Ti3(u)|0) 


V 


(/3(u,0))^|O) 0 

0 (/3(u,0))^|O) 


(36) 


The energy spectrum of the Perk-Schultz spin chain (@ with the (-) grading is obtained using Eq. ( |M| ) and the 

eigenvalues of the transfer matrix (see Appendix [A|) : 


-{v) = 


sin (7 — v) 
sin 7 


^ ^ sin (7 - + u) /sinu 

sin(uj^^ — u) \sin 7 


n 


^ TT sin (7 - vf''^ + v) 

sin(u-uf^) sin(uf^-u) 


f sin 


^,^:^sin(7-u + uf) 

11 ..(2)n 


\sin 7 / sin('(; — ) 


(37) 


with satisfying the BAEs (20), and the fact that the contribution of the chemical potential 

component is given by M{hi — /i2) + Mi{h2 — 7,3] + Eq | 5 ^ . 


B. ABA for the quantum transfer matrix 


The quantum transfer matrix |65L I 68 j is an algebraic object which plays a signihcant role in the investigation of 
finite temperature properties of integrable systems. The importance of the QTM resides in the fact that not only the 
free energy can be derived from the largest eigenvalue but also that different correlation lengths can be characterized 
as ratios of the leading eigenvalues. Pedagogical introductions in the subject can be found in 
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as 


In order to introduce the QTM we need to define two types of L-operators —u), Lj{u,v) € End 

Lj{v,-u)= • (^ 8 ) 


a,6,ai ,fei = l 


a,&,ai ,61 = 1 


where u = — Jsin(ei 7 )^, with N G 4N the Trotter number, and is a canonical basis of operators acting on 
((j^3)(gi(iv+i)^ In the auxiliary space the L-operators can be represented as 


lj{v,-u) = 

^ ai(u, -u)e^^} + /3{v, -u) 
'y+iv,-u)e[^^ 
j+{v,-u)e 


JJ) , AJ) 

^22 ^ ^33 


U) 

13 


a2iv,-u)e^f^ -L/3(u,-u) 

7+(z;, -u)e^'’3^ 


,U) I pU) 
-11 ^ ^33 


7_(u,-u)e^'^2' 


a-i(y,-u)e^^^ + I5{v,-u) 


„0') 

"31 


\ 

,U) 



"32 



u) 

'U) , pU) 
^11 ^ ^22 

/ 


(39) 


Ljiu, v) = 

^ ai(u, v)e^^}+P{u, 

7_(u, 

^ 'y-iu, 


7 +(u, u)e 

a2{u, u)e22^-I-/3 (m, v) 


7_(u, u)e^^2 


e«> 


7+(u, 

OeSi' 

v) 

Jj) , Jj) 

^11 ' ^33 

7+(u, 

1^)4 

) 

jj) 

^32 


asiu, v)e^i^ +I3{u 

, v) 

So) I .0') 

^11 ^ ^22 


\ 

/ 

(40) 


with a canonical basis of operators acting on The monodromy matrix of the QTM, which provides 


(41) 


another representation of the Yang-Baxter algebra ( |32[ ), is given by 

TQ^“(^) = LAr(u, -u)Ln-i{u, v)- - ■ l 2 {v, -u)Li{u, v) . 


Using the representations (391 and (40) we can see that 

/ 


|U) = 


1 
0 

VO 


/o 

0 

yi 


(42) 


N factors 


is a pseudovacuum and 


/ T(3™(u)|U) T^i™iv)\n) 

T^^^iv)\n) = T2^3™(u)|U) r2^™(z;)|U) T2 ^™(u)|U) 

Vt 3 «™(u)|U) Tg™{v)\n} Tg™{v)p 

^ (ai(-y,-M)/3(u,u))^/^|U) {v)\n) 

= 0 (/3(z;,-m)/3(u,u))^/2|U) 

\ 0 0 


Tg™{v)\n) 

Tg™{v)\n) 

(/3(u,-M)a3('u,u))^/2|U) 


The presence of the chemical potential term in the Hamiltonian (18) is taken into account via the transformation 


Q Q 

I 0 0 

0 0 

The QTM is defined as the trace in the auxiliary space of the monodromy matrix 

x9'^^{v) = troT*^™)^) = Ti^i™(u) -L + T§^^{v). 


(43) 

























17 


The existence of the pseudovacuum and the Yang-Baxter algebra ( |3^ ensures that the eigenvalues of for the 

grading (-) can be obtained using ABA with the result (see Appendix [A|) : 


where 


AQ™(^) = Ai(n) + A 2 (n) + A 3 ('!;), 


Ai(u) = -v-u) sin(n - u) ^ g sin (7 - + v) 


Sin 7 sin 7 


\ sin 7 sin 7 J 


A A . / ( 1 ) \ 

sin(v] — V) 

^/2 M _ , , ,,( 1 )^ m (2) 


sin(7 - + ^7) 1 sin(7 - vl ’ + v) 

(l)x 11 .„/,( 2 ) 


sin{v — v) ) sin(u- — I') 


A 3 (u) = sin(7 - m + ^) A sin(7 - u + ) 


sin 7 sin 7 


J.J. . / (2). 

j=i sm[v — Vj j 

with {u7}7i’ solutions of the BAEs Ai(?; 7 )/A 2 (f?^) = -1, = -!)• 


(44) 


VI. FREE ENERGY OF THE (-) PERK-SCHULTZ SPIN CHAIN 

The free energy per lattice site of the Perk-Schultz spin chain is obtained from the largest eigenvalue of the QTM, 
denoted by Ao(u), via the relation /(/ii, /12, /13, / 3 ) = — In Ao( 0 )// 3 . The largest eigenvalue lies in the (A/2, A/2)-sector 
which means that M,Mi = N/2 in Eqs. (44). For our purposes it is useful to change the spectral parameter v —)■ iv. 


Then, the expression for the largest eigenvalue can be written as {N € 4N) 


Ao(i;) = Ai(i;) + A2(u) + A3(n), Xj{v) = (/_(i;)(/+(i;) 


qj-dv + il) qjjv - il) ^ 

qj-id) q,{v) 


,/ 3 /ij 


(45) 


where (/ii,/ 12 ,/ 13 ) = (/ 13 ,/ii,/ 12 ) and 


, . . / sinh(i; ± iu) \ 

(l>±{v) = - - - - 

\ smy 


(/-id 

0 

II 

n^i sinh(z; -vd) 

J = 1,2 

(/+{d 

II 

CO 


(46) 


Using these notations the Bethe equations can be written as \j{vd'^)/\j+i{v)?’) = —1, r = 1, • • • , N/2 . 


„U)\ 


A. Nonlinear integral equations for the auxiliary functions 

We want to obtain an integral expression for the largest eigenvalue of the QTM which can be easily implemented 
numerically. In this paper we are going to use a method which can be understood as the multicomponent generalization 
of the technique presented in |84j (other thermodynamic descriptions of the Perk-Schultz spin chain can be found 
in [87H90) 1. We will introduce a set of two auxiliary functions for which the position of zeroes and poles is known, 
allowing the derivation of NLIEs satisfied by these functions using the Cauchy theorem. In the last step we are going 
to obtain an integral expression for the largest eigenvalue in terms of the auxiliary functions. 

Define 


ai(n) = + d) qid - il) q 2 {v) ^p(hi-h^) 

X2{v) (/-d) qi{v + id q2{v - id 

adv) = gi(^) idv + d) ^0(h2-h,) 

X 2 (v) <p+ (v) qi (v + id q 2 {v - id 


(47a) 

(47b) 


with aj{v) periodic of period in. Each of the equations aj{v) = — 1 has 3A/2 solutions, of which A/2 are the Bethe 

roots and A solutions which are called holes and will be denoted by {vr^'^}di- The distribution of Bethe 

roots and holes characterizing the largest eigenvalue of the QTM for 7 S (0, n/2) is presented in Fig. Let C be a 
rectangular contour with the upper (lower) edges parallel with the real axis through ±^(7 — e)/ 2 , with e —> 0 (please 
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FIG. 8 . Distribution of Bethe roots (‘i ®) a^nd holes (°, □) characterizing the largest eigenvalue of 

the QTM (7 £ (0, 7 r/ 2 )). The contour C contains all the Bethe roots and the poles of the auxiliary functions aj{v) at ±iu. The 
lower edge of the contour C' (see Appendix]^ and next section) coincides with the upper edge of C but has opposite orientation. 


note that the contour is independent of N, the Trotter number). Inside the contour C the function l + Oi(u) (l + 02 (u)) 
has N /2 zeroes at the Bethe roots and a pole of order N /2 at iu{—iu). Then, for v outside of 

the contour, and j = {I, 2} we can introduce 


fAv) = - f — (lnsinh(u — w)) ln(l + aAw))^^ =- f lnsinh(u —in) 

27rl 7c dv 27rz 7c 


a'jiw) 

l + ajiw) 


dw. 


(48) 


where the r.h.s of (48) was obtained by taking into consideration that ln(I + aj{v)) has no winding number due to 
the fact that the number of zeroes inside the contour is equal to the order of the poles at zLiu. We can evaluate fj{v) 
making use of the following theorem: 


Theorem VI. 1 . ^ 91 f Let g{v) be an analytic function inside and on an arbitrary contour C in the complex plane. Let 
4 >{v) be another function which is meromorphic inside and on the contour. Denoting the zeros of 4 >{v) in the interior 
of C by ai, 02, - ■ ■ (with multiplicities ri, r2, • • • ) and the poles by bi,b2, ■ ‘ ‘ (with multiplicities si, S2, • • • ), then 


1 

2TTi 


: 9{v) 


^ r,g{ai) - ^ s^g{hi) , 

ie zeros impales 


with the result 


N 

fi{v) = ln( 3 'i(u) - \n(j)_{v) - — Insiny, 

N 

f2{v) = \nq2{v) - \n(j)+{v) - — Insiny. 


(49a) 

(49b) 


Taking the logarithm of (47) and using (49) we hnd 


Inai(u) = / 3 (/i 3 - hi) + In 


/ (/)-(u - ij) 7>+(u) 
i-f) 4 )-{v) 


+ fiiy- i'y) - fi{v + ij) + f2{v) - f2(v - 17) 


f 2 {v - 17 ) , 

(50a) 

f 2 {v - ij ), 

(50b) 


Eqs. (50) are nonlinear integral equations of convolution type for ai(u) and a 2 (u) which are valid for all N. Using the 
identity 


lim In 
N—yoc 


((+{y) 

((-(v) 


= ij 3 J sinycothu, 
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we can perform the Trotter limit, iV —> oo, obtaining 

J sinh^(* 7 ) 


lnai(z;) = /3(/i3 - hi) - j3 


lna2(?;) = j3{h2 - hi) - /3 — 


sinh V sinh(r> — 17 ) 
J sinh^(* 7 ) 


sinh V sinh(z; + ^ 7 ) 


+ 


+ 


Kq(v — w) ln(l + ai{w)) dw — 

Ki(v — w) ln(l + ai(w)) dw — 


K 2 (v — w) ln(l + 02 (w)) dw, 

(51a) 

Ko{v — w) ln(l + 02 (w)) dw, 

(51b) 


where 


Koiv) = 


1 sinh(2i7) 

27ri sinh(n + ij) sinh(z; — 17 ) 


Ki(v) = 


1 


sinh(i7) 


27ri sinh(t; + 17 ) sinh v 


, K2{v) = 


1 


sinh(z7) 


27ri sinh(ti — 17 ) sinhn 


(52) 


Even though Eqs. ( [M| ) were derived assuming 7 S (0, 7 r/ 2 ) they remain valid also for 7 s ( 7 r/ 2 , tt) by replacing C with 
a similar rectangular contour with the upper (lower) edges situated at ±z( 7 r — 7 — e)/2 with e —?► 0. An observation 
is in order. The above results were derived assuming that v and ± 17 were situated outside the contour C. For v 
on the real axis the integrands that appear in the definition of the functions fi{v ) present an additional pole at v 
which means that in this case we have to add ln(l + 02 (n)) on the r.h.s of Eq. (51a) and ln(l + ai(u)) on the r.h.s of 
Eq. ((sTbl). 


B. Integral expression for the largest eigenvalue 


The free energy of the spin chain is obtained from the largest eigenvalue of the QTM evaluated at zero. Taking into 
account that Ao(n) is analytical in a strip around the real axis then is sufficient to obtain an integral representation 
for Ao (no) with vq inside the strip and then take the limit vq 0. Choosing vq = —iu, for which ^+(rio) = 0, and 
using (Bl) we find 


Ao('yo) = c 


<(>-(^o)g2^^(t^o) 

qi{vo) 


(53) 


Taking the logarithm of (53) and using (B5) we obtain (all constants are denoted by c) 

InAo('yo) = -ln9i(no) + ln(()_('Uo)) + ln9^'*^(no) + c, 

= ln9j'*^(wo) + lng^'‘^(no) - lngi(r>o + * 7 ) - In 92(^0 - * 7 ) - ln(l + ai(no)) + c. 


(54) 


Now we need integral representations for ln( 7 |^^^(tio) and In (rip) with vq close to the real axis. For v inside the 
contour C consider 


2S 


l + ai(it;) zTri 


f ^ J 

- / div — w) -7 —vdw , 

iJc 'l + aiiw) 


(55) 


where on the r.h.s we have used (B4). The latter expression can be evaluated with the use of Theorem VI.1 and 
using that 1 + ai(t!) is periodic of period in and inside the contour C (see Appendix has N zeros which are holes 
and N poles situated at — *7}2i and + 17 }^^ (some modulo in). We find 




<(W) 

1 + ai(ui) 


>3= 

N 


N/2 


N/2 


dw = — ^ d{v — Vj) + ^ d{v — Vj + ij) + ^ d{v — Vj — i"/). 


(56) 


1=1 1=1 
For V inside the contour C, v + ij is situated outside the contour. Then 

N/2 


1=1 


-— f d(v — w + ij)- — ^^ , dw = d(v — + ij) — Nd(v — iu + ij)/2 , 

2ni Jc l + ai(r(;) ^ 


( 57 ) 
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where we have used Theorem VI.l and the fact that inside the contour the function 1 + Oi(w) has N/2 zeros at the 
Bethe roots and a pole of order N/2 at iu. Taking the difference of (561 and (57), integrating by parts w.r.t 

w and then integrating everything w.r.t. v we obtain 

-—: [ [d{v — w) — d{v — w + fy)] ln(l + ai('u;)) dw = lnq 2 {v — fy) — ln( 7 j^^(u) + \n(j)-{v + i'y) + c. (58) 

27ri Jq 

In an analogous fashion we find 

-—: [ [d{v — w) — d{v — w — fy)] ln(l + a 2 ('ic)) dw = lnqi(u + iy) — \nq 2 ^\v) + ln())+(u — iy) + c. (59) 
2Tn Jq 

Replacing hiq[^\vo) and Inq^^^vo) appearing in the second identity of (54) with the expressions which can be derived 
from (1^) and (1^ we find 


In Ao(uo) = ln(())+(uo - i^)4>-{vo + *y))- J Ki{vo - w) ln(l + ai('u;)) dw 

+ J it' 2 (^'o — w) ln(l + a 2 (w)) dw — ln(l + ai(wo)) + c. (60) 


Eq. (60) is in fact valid for vg in a narrow strip around the real axis as a consequence of the analyticity of the largest 

;gration r( 

Ao(u) = 


eigenvalue. Obtaining the constant of integration requires the introduction of 

Ao(u) 


(f>+(v — ij)(/-(v + iy) 


(61) 


which has the property of having constant asymptotics at infinity lim^_,.oo Ao(w) = All the previous 

considerations in this section are also valid for Agiy) which means that 


lnAo(u)=c— / Ki{v — w)\n{l + ai{w)) dw + / Ar 2 ('c — w) ln(l + a 2 (w)) dw — ln(l + ai(w)) 


(62) 


with c a constant and u in a strip around the real axis. Taking the limit u —>■ oo and using lim„_>oo cii(u) = 
and lim„_,.oo a 2 (^’) = we find that 


c = In 


g/3/ii _|_ g/3/i2 _|_ gfihs 
1 _|_ g/3(?i2-/»l) 


= fdhi + In 


1 _|_ gd{h2-hi) _|_ g/3(h3-/ii) 
1 _|_ Q^ih2-hi^ 


(63) 


The last expression for the constant shows that in the scaling limit where /3 = P/S'^ with S O(e^) and 
hi —>■ O(e^), /i 2.3 —>■ O(e^) the term containing the square parenthesis vanishes. Finally, from (61) and using 
lim„_>.o,Ar->.oo ln(^+('y — i'y)4’- {v + *y)) = — d/3 cosy we obtain the integral expression for the largest eigenvalue 

In Ao(0) = c — J/3cosy — y Ar 2 (w) ln(l + Oi(w)) dw + J Ari(w) ln(l + a 2 (w)) dw — ln(l + ai(0)), (64) 


with c defined in (63). This expression, which was derived assuming y S (0,7r/2), remains valid also for y S (7r/2,7r) 
by replacing C with a similar rectangular contour with the upper (lower) edges situated at ±i(7r — y — e)/2. 


VII. CONTINUUM LIMIT 


Having derived an integral expression for the free energy of the spin chain all that remains in order to obtain the 
thermodynamic description of the 2CBG is to perform the continuum limit presented in Section [TV] In the continuum 
limit we can see from Table |T] that y = tt — £ which means that the contour C appearing in Eqs. (51) and (64) has 
the upper (lower) edge parallel with the real axis situated at ±z(7r — y — e)/2. We will denote the upper edge of the 
contour situated at /(tt — y)/2 (the e term is irrelevant for the following discussion) by C+ and by C - the lower edge 
situated at — /(tt —y)/2. For u G C_ , u = x — /(tt —y)/2 with x real, the driving term in the r.h.s. of (51a) is negative 
and equal with 


P{hg-hi)-P 


J " ^ 

J Sin 7 


cosh(x + ^ 7 / 2 ) cosh(x — i^/2) ’ 


/is < /ii, </ > 0, 
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which means that in the scaling limit which implies /3 —> oo , the auxiliary function Ui (u) is very small on C_ and it 
can be neglected (for a rigorous justi ficati on see [92]). In a similar fashion, for u G C+ , u = x + i( 7 r — 7)/2 with x 
real, the driving term in the r.h.s. of (51b) is negative and equal with 


P{h2-h,)-P 


J sin^ 7 


cosh(x + 17 / 2 ) cosh(x — 17 / 2 ) ’ 


/12 < /ii, > 0 , 


which means that the contribution of the auxiliary function 02 (u) on C+ is also negligible. Therefore, in the scaling 
limit Eqs. (51) take the form 

Inai(u) =/ 3 (/i 3 —/ii) — . ' + / Kq{v — w)ln{l + ai{w)) dw — [ K2{v — w) ln{l + a2{w)) dw, 


lna2(u) = /3(/i2 - hi) - (3— 


sinh V sinh(u — 17 ) 
J sinh^(* 7 ) 


(65a) 


sinh V sinh(u + i 


— + / iGi('u — rc) ln(l + ai(w)) dw — / i4ro(^ ~ lii(l + a 2 (w)) dw. 

n) Jc+ Jc- 

(65b) 


We can shift the free argument v and variable of integration to the line +i^l2 for the function ai(u) and to the line 
12 for the function a 2 (u) without crossing any poles of the driving terms obtaining 


In ai(u + 17 / 2 ) = / 3 (/i 3 -hi) - 13 


J sinh^(i 7 ) 


— / Kq(v — w)hi{l + ai{w + i^/2)) dw 


sinh(u + 17 / 2 ) sinh(u — * 7 / 2 ) 

— / K 2 {v — w + — is) \n{l + a 2 {w — i^/2)) dw, 

Jr 


( 66 a) 


In a 2 (u - 17 / 2 ) = / 3 (/i 2 - hi) - (3 


J sinh^(* 7 ) 


— I iGi(u — w — 17 + ie) ln(l + ai(w + 17/2)) dw 

( 66 b) 


sinh(u + 17 / 2 ) sinh(u — i"f/2) 

— / Ko{v — w)lii{l + a 2 {w — i"f/2)) dw. 

Jr 


where we took into account the negative orientation of (7+. Performing similar transformations in the integral 
expression for the largest eigenvalue (64) we find 

In Ao(0) = c—J/? cos 7 + / Ar 2 ('*i'+* 7 / 2 ) ln(l + ai(w+f 7 / 2 )) dw+ / dfi(w— 17 / 2 ) ln(l+ 02 (w— 17 / 2 )) dw—ln(l + ai(0)). 

Jm. Jr 

(67) 

In the continuum limit we have 7 = tt — e, u —>■ Sk/e, w —>■ Sk'/e (we do not need the i factor because we have 
already considered the largest eigenvalue at iv, see the remark before (45)) and 


Ko{v) ^ 


(5 27r 


2 c el 

, Ki{v-i-/)^-- 


S 2tt k{k + ic) 


7^2 (u + *7) —t — 7 


S 27r k(k — ic) 


( 68 ) 


Introducing ai{k) — ai(d/c/e + *7/2), 02 )^) = a 2 (dfc/e — 17/2) and taking the scaling limit in Eqs. ( 66 ) we obtain the 
NLIEs for the two-component Bose gas, Eqs. ([^. 

The grandcanonical potential per unit length for the Bose gas is, see Eq. (24), cj){^,H,i3) = (/(hi, h 2 , ^ 3 ,/3) — 
Eq/L)/S^ with f{hi,h 2 ,h^, (3) = — In Ao(0)/,5 the free energy of the spin chain per lattice site and Eq/L = J cosy—hi 
the zero point energy. In the continuum limit the real part of the driving term on the r.h.s of Eq. (51a) becomes large 
and negative like 0{l/e^) which means that ln(l - 1 - ai(0))/i5 vanishes in this limit. Noticing that Kiidkjt — 17 / 2 ) = 
K2{5k/e + i'j/2) ~ e/27r with dw = 6dk/e and that the second term in the r.h.s of (63) vanishes in the same limit we 
obtain Eq. (|m for the grandcanonical potential. 


VIII. CONCLUSIONS 

In this article we have investigated the density profiles and local density correlation functions of the inhomogeneous 
2CBG at finite temperature. Our results derived using a new and numerically efficient solution for the thermodynamics 
of the uniform system and the local density approximation predict that the polarization at the center of the trap has 
a significant influence on the local correlator throughout the sample this effect becoming more pronounced at strong¬ 
coupling. Even though our analysis of the density profiles did not find any evidence of phase separation this scenario 
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cannot be fully excluded in the region of low-temperatures and small polarizations where our equations (and also the 
TBA result) become numerically unstable as a result of the first-order phase transition. The method employed in 
the derivation of the thermodynamics of the uniform system can also be used in the case of two-component fermions. 
In this case the relevant lattice model is the (—I- -I-) Perk-Schultz spin chain, and will be the subject of a future 
publication [55] . 
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Appendix A: ABA solution of the generalized q = 3 Perk-Schultz model 


In this appendix we present the solution of the generalized q = 3 Perk-Schultz model. The eigenvalues of the 
transfer matrix (33) and quantum transfer matrix (43) can be obtained as particular cases of this general solution. 


We should mention that the transfer matrix (33) was diagonalized in m and the eigenvalues of the quantum transfer 
matrix (43) were conjectured in [88]. For a proof of the results presented below see [94H98] . 


The generalized q = 3 Perk-Schultz model is the set of all linear representations of the Yang-Baxter algebra 

R(u, ic)[T(u) 0 T(u>)] = [T(r)) 0 T(r(;)]R(u, w), 


(Al) 


with the R-matrix defined in (251 and T (v) the monodromy matrix which acts triangularly on a pseudovacuum (highest 
vector) jn) 


/ A(z;)|F!) Ri(^;)|fl) B2{v)\n) \ ( q^i{v)\n) Ri(u)|0) R2(u)|0) \ 

T{v)\n) = Ci{v)\n) D,,{v)\n) D, 2 ivm = o (P 2 (u)|ii) D, 2 {v)\n) . (A2) 

\c 2 {vm D 2 i{vm D 22 {vm) \ o o My)\^)) 


The functions ipj{v),j = 1,2,3 are called the parameters of the model. The monodromy matrices of the transfer 
matrix (31) and QTM (41) satisfy these requirements with pseudovacua (35) and (42) respectively. The parameters 
of the model are given by 


(pi{v) = {ai{v,0))^ , (p 2 {v) = {/3{v,0))^ , (fisiv) = {/3{v,0))^ , (A3) 

in the transfer matrix case and by 

ipi{v) = -u)/3{u,v))^^‘^ , ip 2 {v) = e'^'‘^(/3(u, -u)/3(m, v))^/^ , ipsiv) = -u)a 3 {u,v))^^‘^ , (A4) 

in the QTM case. The eigenvalues of the generalized model are 


A(u) = </5i(u) J|gi(?;y\u) -t ^2{v)Wg2{v,vf’)'^g2{vf\v) + </53(^^) 53(^^,-yf ), 

j=l i=l j=l j=l 


with satisfying the BAEs {gi{v,w) = ai{v,w)/j3{v,w)) 

, (1)n n , (1) (1)n m 


/ (2)\ n m / (2) (2)\ 

_ TT o - 1 ( 7 ;( 2 ) j^y) n ^ 

/ ( 2 ), 11^2 111 / ( 2 ) ( 2 ), 

) i=i j=i92iv) ') 

3^1 


Z = 1, • • • , m . 


(A5) 


(A6) 

(A7) 


In the QTM case, it is preferable to work with a different pseudovacuum (see [94]) which amounts to cycl ic per¬ 
mutations (/ii,/i 2 , Zis) —>■ (^ 3 ,^ 1 ,^ 2 ), (ei; ^ 2 ; £ 3 ); (£ 3 ,£iiE 2 ) of the chemical potentials and grading in (A4). This 

explains the order of the chemical potentials in Eqs. (44). 
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Appendix B: Proof of some identities 


In this appendix we prove several identities which we use in Section VIB First, we will show that 

^ \ ^ ^ (,,) IT • '(2)^ 

A 2 (u) + Xsiv) = c - '[v) = sinh(u - v> '), 


Qiiv) 


(Bl) 


j=i 


with c a constant and the holes corresponding to the equation 02 (u) = —1. Using (45) we have 

(j)-{v)p2{v) 


X2{v) + A3(u) = 


qi{v) q2{v) ’ 


(B2) 


with P 2 {v) = (j)+{v)qi{v + i'y)q 2 {v — * 7 ) 6 ^^^ + 4’+{v — ij)qi{v)q 2 {v + . The function P 2 {v) has the following 

properties: P 2 {v + in) = {—l)^^^‘^P 2 {v) and lim„_,,ooP 2 ('i')/(sinhu)^'^/^ = const. In addition the equation P 2 {v) = 0 
is equivalent with 02 (u) = —1, which means that the zeros of P 2 {v) are the iV/2 Bethe roots and the N holes. This 
shows that 


N/2 


N 


P 2 {v) = c sinh(u - sinh(z; - = c q 2 {v)q 2 ^\v), 


i=i 


f=i 


which together with (B2| proves (Bl). A similar reasoning can be applied to show 

\ r \ \ \ ^+(^) 9 i ^('^) (h), s TT • -m" '(l)^ 

Ai(u) + X 2 (v) = c --, q) '(v) = 11 sinh(u - v/ '), 


q2(v) 


i=i 


(B3) 


where are the holes corresponding to the equation ai(u) = — 1 . 

Consider again the regime 7 S ( 0 , 7 r/ 2 ) for which the distribution of roots and holes is presented in Fig. The 
second useful result which we will prove is 



w) 


1 + aj{w) 


dw = 0 , 


d{v — w) 


— Insinhfu — w ), 
dv 


(B4) 


where C', see Fig. is a rectangular contour with the lower (upper) edges parallel to the real axis through *(7 — e)/2 
and —*(7 — £)/2 + in with £ —:► 0. The lower edge of the contour C at 1(7 — e)/2 coincides with the upper edge of C 
but has opposite orientation which means that the contribution of the J7r-periodic integrand to the integral is zero. 
Then, relation (B4) is proved by noticing that the contributions of the sides parallel to the imaginary axis are also 
zero as a result of limg^i„_>±oo d{v — w) = =Fl and 


a'j{w) a'jiw) 1 
l + aj{w) aj (w) 1 + aj^{w) ’ 


lim 


lim 


lim 


a'(w) 


SRuj—>-±oo 1 -j- 1 hi) ’ sRju—>±00 1 + 02 ^ (rc) 1 + ' SRu)—>-±oo dj {w) 


= 0 . 


Using (Bl) and (B3) and the definition of the auxiliary functions (47) we can also easily derive the following 
identities 


- ln(()_(u) + Ingi(u) - Ingi)!; + iy) - \n.q 2 {v - iy) + \nq[^\v) - ln(l + ai(u)) + Ci = 0 , (B5) 

— ln(()+(u) + lng 2 (i’) — ln(j'i(u + i'f) — lng 2 (u — ij) + Ing^^^(u) — ln(l + a 2 (u)) + C 2 = 0, (B6) 


with ci ^2 constants and arbitrary v. 
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Supplemental Material for EPAPS 

Thermodynamics, density profiles and correlation functions of the inhomogeneous 

one-dimensional spinor Bose gas 


Appendix 1: Algebraic Bethe ansatz for the generalized g = 3 Perk-Schultz model 


Here we solve the generalized q = 3 Perk-Schultz model using the method developed by Gohmann [SS] for 
the gl{l\2) generalized model. The eigenvalues of the transfer matrix (33) and quantum transfer matrix (43) can be 
obtained as particular cases of this general solution. 

The generalized q = 3 Perk-Schultz model is the set of all linear representations of the Yang-Baxter algebra 

R(w, ri;)[T(r;) (g) T(r(;)] = [T(u) (g) T(w)]R(u, w), (1.1) 


with the R-matrix defined in (25) and T (v) the monodromy matrix which acts triangularly on a pseudovacuum (highest 
vector) |0) 


( 1 . 2 ) 


The functions ipj{v),j = 1,2,3 are called the parameters of the model. The monodromy matrices of the transfer 
matrix (31) and QTM (41) satisfy these requirements with pseudovacua (35) and (42) respectively. The parameters 



^ A(?;)|f2) 

Bi(u)|B) 

B2{v)\n) \ 

^ (/7i(u)|B) 

Bi(u)|B) 

B2{vm \ 

T(^;)|B) = 

C'i(^^)lf^) 

Dii{v)\n) 

Di 2 {v)\il) = 

0 

(P2{v)\ft) 

I1i2(u)|B) 


^C'2(^;)|B) 

D2iiv)\n) 

D22iv)\il) J 

1 0 

0 

V?3(w)IT / 


(1.3) 


of the model are given by 

(fiiiv) = , ip 2 iv) = , ipsiv) = , 

in the transfer matrix case and by 

ipi{v) = {ai{v, -u)l3{u,v))^, ip 2 {v) = , ip 3 {v) = e^^^{/3{v,-u)a3{u,v))^^^ , (1.4) 

in the QTM case. We will obtain the eigenvalues of 

t{v) = troT)?;) = A{v) + Dii{v) + D 22 {v ), (1.5) 

using only the Yang-Baxter algebra and the fact that the monodromy matrix acts triangularly on the pseudovacuum. 
It is also important to notice that we do not require that Di 2 {v)\i 2 ) = 0. 

Before we start we need to obtain some auxiliary results which will play an important role in the following consid¬ 
erations. For Q integer and a set of inhomogeneities we consider the following monodromy matrix 


T(?^,{i^j}) = Lq(z),'cq)Lq_i(z;,'cq_i)---Li(w,ui), 


with the L-operators defined as 


( 1 . 6 ) 

(1.7) 


a,b,ai ,bi — l 

where is t he canonical basis of operators acting on (C^)®^'3-i-i)^ Using the formula ~ 

definition (1.7) we can obtain explicitly the elements of the monodromy matrix 

T(i^, eioig ’ ’ ’ - 

as 

TZYZ = E • • • Kill (u V2)Rrb\ (v, vi) • (1.8) 

CQ, " ,C2 = 1 

Then, the elements of the associated transfer matrix f(u, {uj}) = troT(r!, {uj}) are 

= E RcQbQE:'''Q)RcQ_ibQ:iP>'f'Q-i)-"Rc 262 E’'^ 2 )Ra"b;(t',wi). (i.9) 

a,CQ,--- ,C2 — 1 

The previous formulas can be easil y mo dified to the case of different R-matrices and number of lattice sites Q and, 
will be used extensively, especially (1.9), in the following computations. Now we are ready to solve the generalized 
model. 
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1. Bethe ansatz for the first level 


The Yang-Baxter algebra (1.1) can be rewritten in a more explicit form if we use the following block form for the 
monodromy matrix 


T(u) = ( 

C{v) D{v) 

Introducing the matrix 


B{v) = B^iv)), C{v) = 


Ci{v) 

C2{v) 


D{v)= f . 

\D2i{v) D22{v) j 


(li 


X 


0 0 1 
1 0 0 
0 1 0 


V 


( 1 . 10 ) 


I 2 / 


which has the properties: permutes circularly the 5th, 6th and 7th row when multiplied from the left, permutes 
circularly the 5th, 6th and 7th row when multiplied from the right and = X^X = Ig, we can perform the 

following similarity transformation on the Yang-Baxter algebra 0) 


(YR(u,u>)Y^)[Y(T(u)(g)T(w))Y^] = [Y(T(u) 0 T(u>))Y^](YR(u, , 


( 1 . 11 ) 


which can be written as 

^ ai{v,'w) 0 0 0 

0 7+(u,?«)l2 /3{v,w)l2 0 

0 (3{v,w)l2 'y-{v,w)l2 0 


V 


0 


0 


0 R^^^(u,w) J 


^ f A{v) (Si A{w) A{v) (S B{w) B{v) (S A{w) B{v)(SB{w)'\ 
A{v)(SC{w) A{v) (S D{w) B{v)(SC{w) B{v) (S D{w) 
C{v)(SA{w) C{v)(SB{w) D{v) (S A{w) D{v) (S B{w) 

\ C{v) IS) C{w) C{v)SD{w) D(v)sC{w) D{v) S D{w) J 


! A{w) S A{v) A{w) S B{v) B{w) S A{v) B{w)SB{v)\ / a{v,w) 0 

A{w)sC{v) A{w) S D{v) B{w)sC{v) B{w) S D{v) 0 

C{w)sA{v) C{w)SB{v) D{w) S A{v) D{w) S B{v) 0 

\C{w)sC{v) C{w)SD{v) D{w)sC{v) D{w) S D{v) J \ 0 


0 

7+(u,u;)l2 f3{v,w)l2 
I3{v,w)l2 7_(u,w)l2 

0 0 




R^^^(u,iu) J 


( 1 . 12 ) 


where R^^^(t>, w) is the R-matrix of the XXZ spin-chain (modulo a gauge transformation) 


R^^^(u,w) = 


f a2{v,w) 0 0 0 ^ 

0 7 +(u,w) l3{v,w) 0 

0 j5{v,w) j-{v,w) 0 

y 0 0 0 a2{v,w) J 


(1.13) 


For the application of algebraic Bethe ansatz we need the following commutation relations which can be extracted 
from (1.12) 


A{v)A{w) = A{w)A{v ), 

Ba,{v)Ba^{w) = R^^^ ^^\^^^{v,w)Bb^{w)Bb^{v ), 

A(v)Bb{w) = gi{w,v)Bb(w)A{v) — h-{w,v)Bb{v)A{w ), 

Dab^{v)Bb^{w) = gi(u,ui)R^^^ l\ll{v,w)Bc^{w)Dac^{v) - h+{v,w)Bb^{v)Dab^{w) 


where 


R^^^(u,u>) = 


ai(v,w) 


R^^\v,w), g,{v,w) = 


cxiiv.w) 

I3{v,w) 


h±{v,w) = 


l±{ViW) 

(3{v,w) 


(1.14a) 

(1.14b) 

(1.14c) 

(1.14d) 

(1.15) 


t(u)|4>) = A(u)|$), |$)G'H, 


(1.16) 


Our goal is to solve the eigenvalue problem 
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with t(t)) = A{v) + Dii{v) + D 22 {v) and H the Hilbert space of the model. We are going to look for eigenvectors of 
the form 

2 

ai,--- ,0^1 = 1 


where 0 • • • 0 B{vn'^) = (S>J=i can be viewed as a 2” component row vector with entries acting on H, 

F as a, column vector with 2" components belonging to the Hilbert space "H, and are parameters which will 

satisfy Bethe equations. A consequence of commutation relation (1.14b I is the following: 

Lemma 1.1. Let be the transfer matrix of a system of n lattice sites, inhomogeneities and 


R-matrix R^^\v,w) = PFf^' {v,w). Then 


H), 


B{v[^'') 0 • • • 0 B{v^^) = B{v2^^) 0 B{v^^^) 0 • • • 0 B{v^'') 0 . 


(1.17) 


Proof. Starting with Bb,^{v[^'^)Bb 2 {v^^) ■ ■ ■ Bb„{vn'^) and applying successively (1.14b I we find 

BbAv[^^)BbM^)---BbMn^) 


BaAvi^^)BaAvn , 
BaAvi^ABa.Ul^) , 


Comparing with (1.91 modified for the R-matrix and using R ~ ^C 2 bi<^aai we see that 

BbAv[^^)BbAvi^^) ■ ■ ■ BbMnA = })^a2(4'^)Sa3(4'^) ’ ’ ' Ba„ (4'^ ) 

proving ( |1.17[ ) . □ 

Lemma [l.l| also shows that the cyclic permutation B{vl^'^) —)■ followed by the multiplication with the 

matrix U = t^^^ , leaves the product B{v[^'^) 0 • • • 0 B{vA'^) invariant. 

Following the general strategy of the algebraic Bethe an satz w e are going t o act with A(u) -(- Ilii(u) -I- 022 ( 0 ) on 
the eigenstate |$) and then use the commutation relations ( 1.14c| ) and (1.14d) in order to move this operators to the 
right. Due to the fact that in the r.h.s of the commutation relations we have two terms, moving these operators to 
the right will produce 2" terms which we will divide in two categories: wanted and unwanted. The wanted terms 
are those which contain the product B(v[^^) 0 • • • 0 B(vn'’) and the unwanted ones have one of the replaced 

with B(v). In keeping track of the unwanted terms the symmetry described in Lemma 1 1.1 [ will play an important role. 
Satisfying the eigenvalue equation ( |1.16 ) requires that the unwanted terms should vanish, the condition which will 
produce the Bethe equations. 

Action of A(v) on [{S)J=i The wanted term is produced by using only the first term from the r.h.s of 

(1.14c I in moving A(v) to the right obtaining 


[(^ B(v^^^)]A(v) \ u). 


(1.18) 


1=1 


1=1 


The first unwanted term is obtained by using the second term of the commutation relation (1.14c I in order to move 
past B(v[^'^) and then using the first term n — 1 times with the result 


-[B(v) (g) B(v 2 ^^) (g) • • • (g) JJ" ^ 

1=2 

Making use of the U transformation described in Lemma o we can derive a general formula for the fc-th unwanted 
term 

n 

- [B{v) (g) B{vI}IA «)•■•«> B(vAA (g) B{v[^^) g) • ■ • (g) B{vi^}A]U'"~^Aivi^'‘)h-{vl}\v) ]J 5i(i^f \4^^) ■ (1-19) 

1=1 

1/fc 
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Adding (1.181 and (1.19) and introducing the notation 

n 

= B(r;)®B(t>Wi)<g)...0B(^;W)<g)S(z;W)®---(g)S(i;Wi), fc= I,-- - ,n (1.20) 

1=1 

where B{v''p)] = B{v)® B{v^'^) ® ■ ■ ■ ® B{vn^) then 

n n n n n n 


1=1 


1=1 


1=1 


k=l 


1 = 1 


1=1 

iT^fe 


( 1 - 21 ) 

Action of Dii{v) + 1122(11) on [07=i B{v^^'^)]. We will start with Dii{v). The wanted term is obtained by using 
only the first term in the r.h.s of ( 1.14d| ) when moving to the right with the result 

Dn{v)Bt„{v[^^)BtM"^) ■ ■ ■ BbAvi^^) 

n 

1=1 

n 

1=1 


= n 5l(liwf )ila,(lli'^) • • • Sa JllW)Ac(ll)T(°) 

1=1 


( 1 . 22 ) 


The last line of (1.22) was derived using (1.8) with T^°^(u, {u]^^}) the monodromy matrix of a system with n lattice 
sites R-matrix (u, w) and inhomogeneites Introducing 

Zl(o)(u) =il(u)®lf", and T(^)(u, {uf ^}) = H(°)(z;)T(°)(u, {uf ^}), (1.23) 

then T*ic(ii)T^°^ ^^0.’.’.’j)(u, {uj^^}) = T^\^(u, {uj^^}) where we have denoted by {uj^^}) the elements of the 

monodromy matrix T^^^(u, {uj^^}) in the auxiliary space. The computations for 022 ( 0 ) are similar but in this case 
the result involves T22^(i;, {u^^^}). Therefore the wanted term is 

n n n 

(Dniv) + D22{v))[(^B{vl^^)] = [(g) i?(i;f ))] (tW(u, {uf}) + {u f>))) n gi{v,v^^^). (1.24) 


1=1 


1=1 


1=1 


The first unwanted term is obtained using once the second term on the r.h.s of (1.14d) and n — 1 times the first term 
with the result 

D,,{v)B,,iv<i^^)B,M'^)---B,Jvi^^) 

n 

= -h+(u,uW) ■ ■ ■ R^'^ lTbMl\v''2'^)BaAv^l^) ' ' ' B^M^^laMv) 

1=2 
n 

= -h+(v, vZ ) n 91 (y ['^, ) R^'^ Tx ) ■ • • R^'^ S“4 > 4'^ ) R^'^ HZ > 4'^ )Sa, ) • • • i?a„ (vZ^ )Sla, Dr, (v) 

1=2 


IL 

= -h+{v,vHHU9iZi'\yj'^)BaAvHZ---BaAvZHDXv)Z°^rCZl(v[^ 

1=2 


(1.25) 


where the last line was obtained using (1.8) and R^^^ (u^^\ u^^^) = 5c2f)i^aai- Using the definition of the monodromy 

matrix (1.23) the first unwanted term for Zlii(u) is 

n n 

-5°[(g B(vl^Z]Aiiyi\{yj^^})K(y, H 9i(vH\oj^H 

1=1 1=2 
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and a similar expression involving for D 22 {v). The general form of the unwanted term is 

n n 

- }) + {r^ri)) ^+(^^,4'^) n 5 i(4'\4'^) (1-26) 


j=i 


Finally, collecting (1.21), (1.24) and (1.26) we find 


,(l)^ 


1=1 


k=l 1=1 



II 

we find 


n 

[(8)i?(4'4 

/ n n 

i=i 

Vi=i i=i 

71 

71 

4^ ) n > 4^^ 1 

i=i 

i=i 

3¥=k 


2. Bethe ansatz for the second level 


If we introduce the ’’vacuum subspace” [Ml IMl [TOO] . 'Ho & B characterized by the conditions 

A(i;)|4>) = (^i(i;)|$) , C'(i;)|$)=0, |$)G'Ho, 

then we can prove the following 

Lemma 1.2. 'Ho is invariant under the action of D{v). 

Proof. The following relations 

7 +(n, w)B{v) 0 C{w) + /3{v, w)D{w) 0 A{w) = l3{v, w)A{w) 0 D{v) + 7 _(n, w)B{w) 0 C{v ), 

w)D{v) 0 C{w) = fi{v, w)C{w) 0 D{v) + ^-{v, w)D{w) ® C{v ), 

can be obtained from the block form of the Yang-Baxter algebra ( |1.12 ). Using these relations we have 

A{w) ® Il(z;)|<i)) = 73i(w)Il(i;)|$), C{w) ® Zl(w)|<i)) = 0 , 
for all 1$) G T-Lo finishing the proof. 


(1.28) 


(1.29) 

(1.30) 


□ 


A corollary of this lemma is that the linear space spanned by all linear combination of vectors of the form 
• Di 2 {v^m)\ ^) is a linear subspace of 'Ho- 
Using the definition (1.23) and the notation 


T(°Vr; { 4 ^ 4 ) { 4 ^ 4 ) \ 


the elements of T^^^(n, {r'4}) in the auxiliary space are 


(■)(„ („(«()) (131) 

‘'’■U »-(cmK{7')) D<7..{U’})j ' 

_ / A(°)(?;,{4^4)^ii('f')+ C'^°Hr'>{4^4)-Di2(w) B^°'>{v,{v^^^})Dii{v) + D'^°'>{v,{v^^^})Di2{v) \ 

I A(°)(?;,{4^4)^i2(w)+ C'(°)(r'>{4^4)-D22('y) B^°'>{v,{vf'’})Di2{v) + D'^°'>{v,{v^^''})D22iv) J 


where we have used [T*o6('c), {^4J")] “ element (4,4) of (1.12)) and T^°^(u,{u4}) rire both 


representations of the Yang-Baxter algebra with the same R-matrix R*'^^ which means that 

R^^^ (n, w) (v) (g) T^^^ (w)] = (w) 0 T^^^ (u)]R^^^ (u, w). 


(1.32) 
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It follows that trT^^^ (v, can be diagonalized using ABA if we can find a pseudovacuum on which {v, 

acts triangularly. acts on (C^)®" 0 "H and, as we will show below, 


\n) = (g) \n) = 


0 |fi), 


(1.33) 


is an appropriate pseudovacuum. The monodromy matrix (u, is the ordered product of the following 

L-operators 


a,b,ai ,6i = 1 

which can be represented in the auxiliary space as 


Lj{v,vf^) G End , (1.34) 




0) 1. W) U) 

7-(yg^ O') 


(j) 

(1)n 


9i {.v,v^ jeii + „^(^_„(i))e22 


(1.35) 


with a canonical basis of operators acting on (C^)®". Using the representation (1.35) we have 


T 


UU9T\v,v\^^)\n(°^)) 


(1.36) 


and as a consequence of (1.31) 






(1.37) 


which proves that |11) is a proper pseudovacuum. The diagonalization of trT^^^(z;, = A(^^(z;, + 

{uj-^^}) is very similar to the case of an inhomogeneous XXZ spin chain [52]. The following commutation 
relations can be extracted from the Yang-Baxter algebra (1.32) 

A(i)(z;)A(i)(u;) = A(i)(u;)A(i)(u) , 


A''^\v)B^^\w) = g2{w,v)B^^\w)A^^\v) — h-(w,v)B^^\v)A^^\w) , 
D^^\v)B^^\w) = g3{v,w)B^^\w)D^^\v) — h+(v,w)B^^\v)D^^\w) . 


(1.38a) 

(1.38b) 

(1.38c) 

(1.38d) 


We are looking for eigenvectors of the form = B^^'>(v^'^) ■ ■ ■ solving the eigenvalue equation 

trT(i)(u, = A(i)(u)|$(^)) . (1.39) 


Moving A^^)(u) and D^^'>{v) to the right with the help of the commutation relations (1.38c) and (1.38d) we obtain 


(A(i)(u) + Zl(^)('u))[J^ I J^5f2(^^f\'c)A(^^(u) + J^53(w,wj^^)i:)(^^(z;) 


i=i 


i=i 


Vi=i 


1=1 


/ 


\ 


- XI n )] ^'^)W92{vf’ , Wfc ) [ v ]^’) + h+{v, v)^’ ) 53 (4^^ > (w^) 


k=l 


1=1 


V 


1=1 

l#fc 


1=1 

l#fc 


/ 


(1.40) 
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Using and HLi 5i then 

O'! (17,1^ - ) J J 


TU 0^2 («,u 
h~l cKi 

is an eigenvector of trTi^i(r;, with eigenvalue 


/ (i)\ 

a2(v,v) ') 




(1.41) 


(1.42) 


A(i)(t>)=y 2 (t’)n ’ m. +‘^3(^)n^i)n53(^> 

j=iai{v,v) ') J =1 i=i 

if the first set of Bethe equations is satisfied. 

/ (2)\ n m / (2) (2)\ 

*^2(4 ) „TT - 1 / (2) k-1--- m 

. ( 2 )^ 11^2 .I’i U 1 . ( 2 ) ( 2 ). ’ * 

/ 2=1 j—1 / 

3^k 

As a consequence of the invariance of the vacuum subspace Hq under the action of D{v), i3(^i(up^) • • • (vm^ )\0,) 

is a column vector with 2” rows having vectors in Ho G H. Therefore, i® an eigenvector of 

A{v) seen as an operator on (C^)®" 0 H 

A(u)S(i)(u4) • • • B(^\v^^^)\n) = <^i(u)i?(i)(uf) • • • B<^^\vl^'>)\n) . 

An immediate consequence of this relation and ( |1.27 l, (1.39) is that 

n 

[(g) B{vf^)]B^^\vf^) ■ ■ ■ = B,,(4'^) • •. G H 

i=i 

is an eigenvector of A(u) with eigenvalue 

n 71 m m 

A{v) = </ 5 i(w) +7^2(^^)n52(^:'*^4) n52(4^\w) +(/j 3 (z;) ] 4 g 3 (u,wf^), 


(1.43) 


i=i 


i=i 


i=i 


if the second set of Bethe equations 


/ (1), n , (1) (1), m 

T’lW ) _ rr 52(u ) TT ^ ,(2) „ ,(1) 


n 


7>2(v'P) i=i 4^^) M 

i^k 




3 ’ 


0 , 1 = 1,- 


(1.44) 


is satisfied. 


3. Choosing a different pseudovacuum 


of the chemical potentials (hi, ^ 2 , / 13 ) 
how we can obtain (44) by employing a different pseudovacuum, 
monodromy matrix (41) on the pseudovacuum 


The formula (|44|) for the eigenvalues of the QTM differs from (1.431 with parameters (1.4) by a circular permutation 

{ho, hi, ^ 2 ) and grading (ei, £ 2 , £ 3 ) —f (£ 3 , £i,£ 2 )- In this section we will show 

Making use of (39) and (40) the action of the 


is given by 


T^^^{v)\n) = 


\n) = 



/T,T“(^;)|f2) 


QTM/ 

Ti™{vm 

Tg™{v)\h) 


?2™iv)\n) 

QTM, 

22 


t: 


T^A^''^{v)\n) 


t: 


QTM 


N factors 


(^)|U) \ 


Tg™iv)\n) 


13 

QTM ( 

23 

Tg™{v)\n)) 




e^hi (^( 1 ;^ r;))^/^ |U) 


Tg™{v)\n) 


(1.45) 


(1.46) 


0 

nQTM, 
■^31 


g/3/»2 {I3{v^ —u)a2{u,v))^^‘^\^) 




T) 


QTM 


32 


iv)\n) 


{ao{v, —u)f5{u, u))-^/^|U) 


























Consider now the similarity transformation of the Yang-Baxter algebra (32) 


(xni2n23R(r;,u;)n^3nr2^^)[^ni2n23(T«™(z^) 0 T«™(u;))n^3n?2Y^] = 


[xni2n23(T^™(r^) 0 " (u;))n^3n^2Y^^ ](Yni2n23R(?;, «;)n^3n^2^'') 


rQTM 


with X defined in (1.10 1 and ni 2 = 7 ri 2 C) 7 ri 2 , 1123 = 7r23 <8 7123 , 


7^12 


^0 1 o\ 
10 0 , 
Vo 0 1 ) 


7r23 


1 0 o\ 
0 0 1 
0 10 / 


(1.47) 


(1.48) 


The matrix 7112 ( 7123 ) permutes the states 1 and 2 (2 and 3) in the auxiliary space and 111211^2 = ni 2 ni 2 = I 9 , 112311^3 = 
n23n23 = I 9 . Then, we have 


7ri27r23T'^™ (^')7^ J37rf21 


/ Tg^^iv)\n) Tg™ivm 

Tfi™{vm T^i™{v)\n) T^™{v)\n) 
\Ti™{v)\n) Ti™{v)\n) Ti™{v)\n} 


with 


/ v?i(77)|f^) T3«™(u)|ii) \ 

0 (^' 2 ( 77 ) 111 ) T^™{v)\n) , 

\ 0 0 i^ 3 (u)|fi) / 


ipi{v) = e^^^{a^(v,-u)l3(u,v))^/'^ , ip 2 {v) = e^^7(/?(u,-M)/?(u,z;))^/^ , 


ip^(v) = e^^^{l3{v,-u)a2(u,v))^/'^ (1.49) 


Together with 


where 


Yni2n23R(77,7c)n^3nf2Y^ 


^ 03 ( 77 , w) 0 0 

0 7_(w,ri7)l2 P{v,w)l2 

0 j3{v,w)l2 7 +( 77 , 7 c)l 2 

V 0 0 0 


0 

0 

R^^Vt’jW) / 


R^^^( 77 , w) 


/ 01 ( 77 , 777 ) 0 0 

0 7 +( 77 , 777 ) l3{v,w) 

0 l3{v,w) j-{v,w) 

V 0 0 0 


0 \ 

0 

0 

02(77,777) / 


(1.50) 


(1.51) 



























